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ABSTRACT 

The acceleration of a spherical dust particle caused by an interstellar gas flow depends 
on the drag coefficient which is, for the given particle and flow of interstellar gas, a 
specific function of the relative speed of the dust particle with respect to the interstellar 
gas. We investigate the motion of a dust particle in the case when the acceleration 
caused by the interstellar gas flow (with the variability of the drag coefficient taken into 
account) represent a small perturbation to the gravity of a central star. We present the 
secular time derivatives of the Keplerian orbital elements of the dust particle under 
the action of the acceleration from the interstellar gas flow, with a linear variability of 
the drag coefficient taken into account, for arbitrary orbit orientation. The semimajor 
axis of the dust particle is a decreasing function of time for an interstellar gas flow 
acceleration with constant drag coefficient and also for such an acceleration with the 
linearly variable drag coefficient. The decrease of the semimajor axis is slower for the 
interstellar gas flow acceleration with the variable drag coefficient. The minimal and 
maximal values of the decrease of the semimajor axis are determined. In the planar 
case, when the interstellar gas flow velocity lies in the orbital plane of the particle, 
the orbit always approaches the position with the maximal value of the transversal 
component of the interstellar gas flow velocity vector measured at perihelion. 

The properties of the orbital evolution derived from the secular time derivatives are 
consistent with numerical integrations of the equation of motion. The main difference 
between the orbital evolutions with constant and variable drag coefficients is in the 
evolution of the semimajor axis. The evolution of the semimajor axis decreases more 
slowly for the variable drag coefficient. This is in agreement with the analytical results. 
If the interstellar gas flow speed is much larger than the speed of the dust particle, 
then the linear approximation of dependence of the drag coefficient on the relative 
speed of the dust particle with respect to the interstellar gas is usable for practically 
arbitrary (no close to zero) values of the molecular speed ratios (Mach numbers). 
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1 INTRODUCTION 

Recent observations of debris disks around stars with asym- 
metric morphology caused by the mo tion of the stars 
through clouds of i nterstellar matter (IHines et al.l 120071 : 
iManess et all 120091 ; iDebes. Weinberger fc Kuchnen 12009? ) 
have presented evidence that the motion of a star with re- 
spect to a cloud of interstellar matter is a common pro- 
cess in galaxies. The orbital evolution of circumstellar dust 
particles is investigated for many decades. From accelera- 
tions caused by non-gravitational effects accelerations due 
to the electromagnetic and corpuscular radiation of the star 
are most often taken into account. They are u sually de- 
scribed by the Poynting-Robertson (PR) effect (|Povntingl 
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1903; Robertson! 1937) and radial stellar wind (|Whipplel 
1959 ; iBurns. Lamv fc Sotdll979l : lGustafsoijli994 ). respec- 



tively. The acceleration acting on a spherical body moving 
through a gas, derived under the assumption that the ra- 
dius of the sphere is small compared with the mean free 
J iath of the gas, was published a re latively long time ago 
Baines. Williams &: Asebiomdll965r ). However, the first at- 
tempt to describe the orbital evolution of circumstellar dust 
particles under the actio n of an interstellar gas flow was 
made relatively recently (Schcrcr 2000). Scherer has calcu- 
lated the secular time derivatives of the particle's angular 
momentum and the Laplace-Runge-Lenz vector caused by 
the interstellar gas flow. When the interstellar gas flow ve- 
locity vector lies in the orbital plane of the particle and the 
particle is under the action of the PR effect, radial stel- 
lar wind and an interstellar gas flow, the motion occurs in 
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a plane. In this planar case the secular time derivatives of 
the semimajor axis, the ecc entricity and argum ent of the 
perihelion were calculated in lKlacka et al.l |2009). The sec- 
ular time derivatives of all Keplerian orbital elements under 
the action of an interstellar gas flow with constant drag co- 
efficient for arbitra ry orbit orientation were calculated in 
(jPastor et al.ll201ll ). In this paper, it is analytically shown 
that the secular semimajor axis of the dust particle under 
the action of an interstellar gas flow with constant drag coef- 
ficient a l ways decreases. This result contradicts the results of 
ISchererl (120001 ). He came to the conclusion that the semima- 
jor ax is of the dust particle increases exponentially ([Schererl 
120001 . p. 334). The de crease of the semimajor a xis was con- 
firmed an alytically by Belvaev fc Rafikov (2010) and numer - 
ically bv iMarzari fc Thebaultl (120111 ) and lMarzaril (J20121 ). 
iBelvaev fc Rafikovl (120101 ) investigated the motion of a dust 
particle in the outer region o f the Solar system behind t he so- 
lar wind termination shock. IBelvaev fe Rafikovl (|2010h used 
an orbit-averaged Hamiltonian approach to solve for the or- 
bital evolution of the dust particle in a Keplerian poten- 
tial subject to an additional constant force. The problem 
which they solved is known in physics as the classical Stark 
problem. If the speed of the interstellar gas flow is much 
greater than the speed of the dust grain in the stationary 
frame associated with the central object, and if the speed of 
the interstellar gas flow is also much greater than the mean 
thermal speed of the gas in the flow, then the problem of 
finding the motion of a dust particle under the action of the 
gravity of the central object and of the interstellar gas flow 
reduces to the classical Stark problem. The secular sol ution 
of Stark problem presented in lBelvaev fc Rafikovl 1)20101 ) was 
confirmed and generalis ed using a different perturbative ap- 
proach in lPastorl (J20121 ). 

In this paper, we take these studies a step further 
by taking into account the variability of the drag coef- 
ficient in the acceleration caused by the interstellar gas 
flow. An explicit form of the dependence of the drag co- 
efficient on the relative speed of the dust particle with re- 
spect to th e interstellar gas was derived already in the paper 
Baincs, Williams & Ascbiomo (1965). Belvaev & Rafikov 
(J2010j) calculated the secular time derivative of the semi- 
major axis using the linear term in the expansion of the 
drag coefficient into a series in the relative speed of the dust 
particle with respect to the interstellar gas. In this paper, we 
calculate the secular time derivatives of all Keplerian orbital 
elements with the linear term in the expansion taken into 
account. We compare the orbital evolution with the con- 
stant, linear, and explicit dependence of the drag coefficient 
on the relative speed of the dust particle with respect to the 
interstellar gas. 



2 SECULAR EVOLUTION 

The acceleration of a spherical dust particle caused 
by a flow of neutral gas can b e given in the form 
(jBaines. Williams fc Asebiomdll965l ) 



dr 



— = — 2jcDi 7^ |i7— {/f| (v — Vf) 



dt 



(1) 



The sum in Eq. (TTJ runs over all particle species i. vf is the 
velocity of the interstellar gas flow in the stationary frame 
associated with the Sun, v is the velocity of the dust grain, 
cm is the drag coefficient, and 7j is the collision parameter. 
The drag coefficient can be calculated from 



CDi(Si) 



0? \Si 2s^ 



(2) 



where erf(s^) is the error function erf(si) = 2/^/n ff* e _t dt, 
Si is the fraction of impinging particles specularly re- 
flected at the surface (for the resting particles, there is 
assumed diffuse reflec tion) (jBaines. Williams fc Asebiomd 
1 196a ; iGustafsonl Il994l ). Td is the temperature of the dust 
grain, and T is the temperature of the ith gas component. 
Si is defined as a molecular speed ratio 



nii 
2kT~ 



U 



(3) 



Here, m,i is the mass of the neutral atom in the ith gas 
component, k is Boltzmann's constant, and U = \v — vf\ 
is the relative speed of the dust particle with respect to 
the gas. The dependence of the drag coefficient on Si for 
specular (Si — 1) and diffuse (Si — 0) reflection is depicted 
in Fig, [l] For diffuse reflection, we assumed that Td = Ti 
(jBaines. Williams fc Asebiomd 1 19651 ). The drag coefficient 
is approximately constant for Si 3> 1. However, if the in- 
equality Si ;§> 1 does not hold and changes of the relative 
speed U during orbit are not negligible, then cm depends on 
U and can not be approximated by a constant value. There- 
fore, in this case is necessary take into account dependence 
of cm on the relative speed U. For the primary population of 
the neutral interstellar hydrogen penetrating i nto the Solar 
syste m we obtain si = 2.6 using Ti = 6100 K dFrisch et alj 
|2009| ) and U = \v F \ = 26.3 km s" 1 (|Lallement et al J 12005ft 
in Eq. @. Because inequality si 3> 1 does not hold for this 
value of the molecular speed ratio (Mach number), variabil- 
ity of the drag coefficient can have interesting consequences 
also in the Solar system. For the collision parameter, we can 
write 



mi 
Ji — ni — A 
m 



(4) 



where n; is the concentration of the interstellar neutral 
atoms of type i, and A = ttR 2 is the geometrical cross sec- 
tion of the spherical dust grain of radius R and mass m. For 
Si C 1, or, more precisely, if si is negligible in comparison 
with Si is possible to show that 
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Figure 1. Dependence of the drag coefficient cpn on the molec- 
ular speed ratio Si for the cases of specular and diffuse reflection 
(see text). 
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©. We have 



CDi{Si) 



where 

SOi = 
COi = 



CDi(SQi) + 



CDi(S0i) + 



dc Di 
dsi 

dc Di 
dsi 



CO, 



2kTi 

CDi(soi) 

f dc Di 
V dsi 



V ■ Vp 
Vf 



Vf , 



(si — SOi) 



mi irr \ 

2m iu ~ VF) 



(10) 



2kTr 



(11) 



We can rewrite Eq. © using these two approximations into 
the following form 



As a consequence acceleration of the dust particle is 
IJBaines. Williams fc Asebiomdli965l ) 



fir 
It 



E 



8 1 
3 0r 



! + (!-&) 



J_d \ 7T 

Ti 



1i 



2kTi , 

x \l (v - vf) ■ 



(6) 



Hence, for s; < 1 acceleration depends linearly on the rel- 
ative velocity vector v — vf- The case Si <^ 1 will be no 
further discussed in this parer. 

We want to find the influence of an interstellar gas flow 
on the secular evolution of a particle's orbit. We assume 
that the dust particle is under the action of the gravitation 
of the Sun and the flow of a neutral gas. Hence, we have the 
equation of motion 



dv 

It 



/'• 



E 



cm Ji \v — v F \ (v — Vf) 



(7) 



where /i = GMq, G is the gravitational constant, Mq is the 
mass of the Sun, r is the position vector of the dust particle 
with respect to the Sun, and r — \r\. 

We will assume that the speed of the interstellar gas 
flow is much greater than the speed of the dust grain in the 
stationary frame associated with the Sun: 



\v\ = v <S \vf\ — vf ■ 
Therefore, we can write 



(8) 



dv 

Tit 



+ 



i 



2 

V F 



v vf 
vf vf 



Co, 



-VF 



V ■ Vp Vp 

v F vf 



(12) 



This equation allows using the perturbation theory of celes- 
tial mechanics to compute the secular evolution of the dust 
particle under the action of the interstellar gas flow. For the 
secular time derivatives of the Keplerian orbital elements 
caused by the interstellar gas flow, we finally obtain (see 
Appendix [All 



da \ s—^ 2 pP 1 1 



I da\ 



X ( 1 H VF 

COi 



(J 2 - S 2 ) 



1 - VI - e 2 



de\ 

dt) 
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COi Ji V F 



(13) 



37 a(I 2 - g 2 )(l - e 2 ) 
2 VFe 3 



Co, 



l + ^v F )(l- C --VT- 



(14) 
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IV — Vf\ 



VF 1 



■!'- 



v ■ Vf 

~vl~ 



2 v ■ vf 



(9) 



In the above equation, we have considered only the terms 
to the first order in v/vf. Using this approximation, we can 
also approximate changes in the drag coefficient cm in Eq. 
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Sesinui a /, ki 

t 2 1 + — v F 

1 — e z vf V c 0i 



x (S cos uj — I sin uj) 



(15) 



d\l\ y^ 

~Ht I ~ 2-^ 



CQi Ji VF C 

2 sini 



3e sin uj 
1-e 2 



— [ 1 H -«p ) (Scoscj — J sin a;) 

v F V c o» 



(16) 









/ di \ _ ^ c 0i 7» «f C 
\ di / ^ 2 

* ' i 


\ll 


3e cos uj 
1-e 2 


a ( ki \ 
+ — 1 + — VF 

VF V C 0i J 


(Ssina; + I cos uj) 


where p = a(l — e 2 ), 




a = . 







VF 

and the quantities 

S = (cosfl cosw — sinfi sin a; cosi) vpx 

+ (sin cos u> + cos Q sin w cos i) vfy 

+ sin a; sini ufz , 
I — (— cosil sina; — sinfi cosw cosi) dfx 

+ (— sinfi sincj + cosfi cos a; cosi) vfy 

+ cos a; sini vfz , 
C — sin 57 sin i vfx — cos fi sin i ufy 

+ cos i vfz , 



, (17) 



(18) 



(19) 



are the values of A — vf ■ e*R, B — vf ■ St, and C = vf ■ ejv, 
at the perihelion of the particle's orbit (/ = 0), respectively. 
For a complete solution of this system of equations for a = 
(constant force) , we refer the reader to IPastorl (|2012i ) . 



3 DISCUSSION 

C — for the special case when the velocity of the interstel- 
lar gas, vf, lies in the orbital plane of the particle. In this 
planar case, we find that the inclination and the longitude 
of the ascending node are constant. 

Eq. (|15[) implies that the argument of the perihelion 
is constant in the planar case (C = 0) and if the orbit's 
orientation is characterised by S — 0. 

The dependence of the drag coefficients on the relative 
speed of the dust particle with respect to the interstellar gas 
is demonstrated by the presence of terms multiplied by ki 
in Eqs. (|13[l - (|17p . It is convenient to define a new function 
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CQi 



(20) 
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Figure 2. Dependence of Qi on SQi for the case of specular and 
diffuse reflection. 



In order to find the influence of a non-constant drag coeffi- 
cients on the evolution of the particle's orbit we will analyse 
the properties of this function. We can write, see Eqs. (|11[) . 



= 14 
1 

CQi 



dc D , 
dsi 



■so; 
COi 



2s 3 

zb 0i 



q2 i 4 „4 , erf(ao*) 



+ 



(21) 



The graph of gt for the case of specular (Si — 1) and diffuse 
(Si — 0, T d — Ti) reflection is depicted in Fig. [5] The func- 
tion gi is an increasing function of SQi for so; G (0, oo) (see 
Appendix[B]). lim Sgi ^o Qi — and linis^-joo gi = 1. Hence, 
we can conclude that gt £ [0, 1]. 

Eq. (|13l) can be rewritten in the following form. 



da \ 
It / 



— } 2 a cot 7< v F a / — a 



1 + 



X ( H V F 

CQi 



1- Vl^ 



(/V 



1 - e 2 + S 2 



(22) 



Thus, the semimajor axis is a decreasing function of 
time. This r e sult, when ki — 0, was already obtained in 
Pasto r et al.l (12011 ), and generalised to the case ki 7^ in 
Belvaev fc Rafikovl ((2010J). If we use the properties of gi, 
then from Eq. (|22[) we can conclude that the dependence of 
the drag coefficients on the relative speed of the dust parti- 
cle has a tendency to reduce the decrease of the semimajor 
axis caused by the interstellar gas flow. 

In order to find the orbit orientation with minimal and 
maximal decrease of the semimajor axis, we will analyse the 
second term in the square brace in Eq. (I22|l . 
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Figure 3. An orbit orientation with minimal decrease of the semi- 
major axis. 



— |H -vf 

C i 



1-vT 



(i 2 VT 



e 2 + S" 



(23) 



Because the terms multiplied by S 2 and I 2 are both 
positive, we obtain a minimal value of <f> when S = and I 
— 0. Therefore, if the orbital plane is perpendicular to the 
interstellar gas flow velocity vector, then the decrease of the 
semimajor axis is minimal (Fig. [3]). From Eq. (|22[1 . we obtain 
for S = and 7 = 0, 



da \ 

~dl) 



E 2 « 



coi 7i vf 



(24) 



The value of the minimal decrease is proportional to the 
semimajor axis and independent of the orbit eccentricity. 

Because the terms multiplied by S 2 and I 2 are both 
positive, we obtain a maximal value of <j> when C = 0. If C 
— 0, then S 2 + I 2 — v F . Using this, the value of <j> can be 
written as 



= t[H -Vf 

COi 



l-VT 



x \v%Vl - e 2 + S 2 (l -y/l- e 2 )l 



(25) 



Here, v F is constant. Therefore, we obtain the maximal 
value of 4> for an orbit orientation characterized by S 2 = 
Vp. Hence, the maximal value of <f> is 



= ( 1 + ^VF 

COi 



1- vT^ 



= 1 



— vf h(e) . (26) 

COi ■ 



Therefore, if the interstellar gas flow velocity vector is paral- 
lel to the line of apsides, then the decrease of the semimajor 
axis is maximal (Fig. |4|). For a given orbit, the maximal 
decrease of the semimajor axis is 



da \ 
~dl) 

I mas 



2J 2 a c i 7 t 
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vf 


1- Vl-e 2 




e 2 





I - i 1 + — VF 



(27) 



The function h(e) defined in Eq. (|26p is an increasing func- 




Figure 4. An orbit orientation with maximal decrease of the 
semimajor axis. 



tion of the eccentricity (see Appendix [Cj . Therefore, the 
decrease of the semimajor axis is maximal for e = 1. 

For the secular time derivatives of S, I, and C, we ob- 
tain from Eq. Jl9j) and Eqs. [fTS]). (fj"6)) and (JTTJ) , 



dS_ 
dt 



COi 7i v F S 



dl\ 

~dl J 



dC 

~~dt 




(28) 



(29) 



(30) 



Eqs. (|28[) - (|30[) are not independent, because S{dS/dt) + 
I{dl/dt) + C{dC/dt) = always holds. Eqs. p8 )l -p0 |l . to- 
gether with Eqs. (I13|l and (|14|l . represent the system of equa- 
tions that determines the evolution of the particle's orbit in 
space with respect to the interstellar gas velocity vector. All 
orbits that are created from rotations of one orbit around 
the line aligned with the interstellar gas velocity vector and 
going through the centre of gravity will undergo the same 
evolution determined by this system of equations. If a is 
small and / and e are not close to zero, we can use the follow- 
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ing a pproximate solution for S, I, and C fsee lPastor et al.l 
[203). 



e 



C; 



V 



VT 



(31) 
(32) 



4 NUMERICAL RESULTS 

4.1 Accelerations influencing the dynamics of 
dust grains inside the heliosphere 

For a correct description of the motion of micron-sized dust 
particles inside the heliosphere, solar electromagnetic radia- 
tion and the solar wind must also be considered. 



and 

171 ft 



LP 

c 2 



V 2 



(33) 



where U and V are some constants. 

Now, we want to find the evolution of the orbit position 
in the planar case. For this purpose we can use Eq. (12911 , 
which determines the time evolution of /. Eq. ((29} implies, 
for the planar case (C = 0), 



\~di 



E 



cot "/i Vf S 



1 + — v F ) 6(e) 

COi 



vf 



Here, 



6(e) 



e 4 -6e 2 +4-4(l-e 2 ) 3 / 2 



(34) 



(35) 



The function 6(e) is a dec reasing function of eccentricity for 
e G (0, 1] fsee lPastor et alj|201ll . Appendix B). The function 
6(e) attains values from lim e _>o 6(e) = — 0.5 to 6(1) = — 1, 
for e G (0, 1]. Since we have assumed that v <C Vf, see Eq. 
JSJ, we have for the maximal speed of the dust particle in 
the perihelion of the particle's orbit, see Eq. (|A4[) . 



(1 + e) < v F 



Hence 



(l + e)< 1 



(36) 



(37) 



For / > 0, one always has (dl/dt) > 0. Therefore, we will 
assume that I < 0. For negative I, we can write 

— ( 1 + —Vf) 6(e) < — 6(e) < -a 6(e) < a < - , (38) 
vf \ Coi J ' vf e 



as 1 + ki vp/coi ^ 1, see the discussion after Eq. (|21[) , — I 
^ Vf, 6(e) G (— 0.5, — 1], and a< 1. If we rearrange (|38l) , 
then we come to the conclusion that (dl/dt) > also for I 
< 0. Therefore, in the planar case, I is always an increasing 
function of time. Thus, in the planar case the orbit rotates 
into position with a maximal value of /. In this position, the 
line of apsides is perpendicular to the interstellar gas flow 
velocity vector. 



4-1.1 Electromagnetic radiation 

The acceleration of a dust particle with spherically 
distributed mass caused by electromagnetic radiation 
to t he first order in v/c is given by the PR ef- 
fect jPovntind Il903l: iRobertsonl Il937l: jWvatt fc Whipple 
ll95d;lBurns. Lamv fc Soterll 19791 ; lKlackal2004l ; lKlacka et al 
l2012bh : 

dt P r 2 



, v ■ eR \ _ v 

1 e R 

c c 



(39) 



where Sr — f/r and c is the speed of light in vacuum. The 
parameter /3 is defined as the ratio of the electromagnetic 
radiation pressure force and the gravitational force between 
the Sun and the particle at rest with respect to the Sun 



P = 



3 Lq Q pr 

16 ty c n R g 



(40) 



Here, Lg is t he solar luminosity, L Q = 3.842 x 10 6 W 
(|Bahcallll2002f ). Q' pr is the dimensionless efficiency factor for 
radiation pressure integrated over the solar spectrum and 
calculated for the radial direction (Q' pr = 1 for a perfectly 
absorbing sphere), and g is the mass density of the particle. 



4-1.2 Radial solar wind 

Acceleration caused by the radial solar wind t o the first or- 
der ot v/c and the first order of v/u is given by IjKlacka et al.l 
2012a, Eq. 37): 



dv r\ u (i 

dt Q' c r 2 



v ■ eR \ „ v 
I I en 



(41) 



Here, u is the speed of the solar wind with respect to the 
Sun, u — 450 km/s. r\ is the ratio of solar wind energy to 
electromagnetic solar energy, both radiated per unit of time 



2 N 

4irr ii \— * 

! ^ ? ZV 

! = 1 



^O 



where m s 



and n. 



(42) 



i — 1 to N, are the masses and 



concentrations of the solar wind p articles at a distanc e r 
from the Sun. r] = 0.38 for the Sun iJKlacka et al.ll2012al '). 



4-1.3 Acceleration caused by solar gravity, solar radiation, 
and interstellar gas flow 

In order to find the acceleration of the dust particle inside 
the heliosphere, we can sum the gravitational acceleration 
from the Sun, the acceleration from the PR effect Eq. (|39|) , 
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the acceleration from the solar wind Eq. (|41|) . and the ac- 
celeration from the interstellar gas Eq. (JT}. 



dv 
~di 



(1 - P) e R 



P^(l+ V 



v ■ e« _ v 

cr + - 

c c 



— ^^CDi-fi \v — V F \ (v — V F ) 
i 

Here, it is assumed that (r}/Q' pr )(u/c) <C 1. 



(43) 



4.2 Comparison of the solution of the equation of 
motion with the solution of the system of 
equations constituted by the secular time 
derivatives of the Keplerian orbital elements 

We want to compare the solution obtained from Eq. ()43p 
with the solution of the system of equations constituted 
by the secular time derivatives of the Keplerian orbital el- 
ements. To do this, we need to add to the right hand sides 
of Eqs. (|13[) - (|17[) also the secular time derivatives of the 
Keplerian orbital elements caused by the PR effect and the 
radial solar wind. Therefore, we solved the following system 
of equations (|Wvatt fc WhippleH T950: Klac ka et alj|2012al ) 
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As the central acceleration, we used the Keplerian acceler- 
ation given by the first term in Eq. (I43[) . namely — fj,(l — 



P)e R /r . This is denoted by the subscript P in Eqs. ((44)) - 
Q48p . In the interstellar gas flow, we have taken into ac- 
count the primary and secondary populations of neutral 
hydrogen atoms and neutral helium atoms. The primary 
population of neutral hydrogen atoms and neutral helium 
atoms represent the original atoms of the interstellar gas 
flow which penetrate into the heliosphere. The secondary 
population of neutral hydrogen atoms are the former protons 
from the interstellar gas flow that acquired electrons from 
interstellar H° between the bow shock and the heliopause 
(|Frisch et all 120091 ; lAlouani-Bibi et al.ll20liT ). We adopted 
the following parameters for these components in the inter- 
stellar gas flow, m = 0.059 cm" 3 and 7 \ = 6100 K for th e 
primary population of neutral hydrogen (|Frisch et al.l l2009), 
ri2 = 0.059 cm - and T2 = 1 6500 K for the sec ondary pop- 
ulation of neutral hydrogen (JFrisch et al.ll2009r i and finally 
TI3 = 0.015 cm" 3 and T3 = 6300 K for the neutral helium 
IjLallement et al.ll2005T ). We have assumed that the interstel- 
lar gas velocity vector is equal for all components and identi- 
cal to the velocity vector of the neutral helium entering the 
Solar system. The neutral helium enter t he Solar system 
with a speed of about v F = 26.3 km s _1 (|Lallement et al.l 
120051 ). and arrive from the direction of A ec ; = 254.7° (he- 
liocentric eclipti c longitude) and p ec i = 5.2° (heliocentric 
ecliptic latitude; lLallement et al.l 12005). Thus, the compo- 
nents of the velocity in the ecliptic coordinates with the 
a;-axis aligned towards the actual equinox are v F — — 26.3 
km/s [cos(254.7°)cos(5.2°), sin(254.7°) cos(5.2°), sin(5.2°)]. 
We want also to demonstrate the influence of a variable drag 
coefficient on the secular orbital evolution of the dust parti- 
cle's orbit. Therefore, we solved Eq. (|43|) and the system of 
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Figure 5. A comparison of the solution of the equation of motion (solid lines) with the solution of the system of differential equations 
constituted by the secular time derivatives of the Keplerian orbital elements (dashed lines). The solutions with variable (black lines) and 
constant (grey lines) drag coefficients are compared. 



Eqs. (I44[) - (|48|) in two cases. One with variable drag coeffi- 
cients and one with constant drag coefficients. The variable 
drag coefficients for Eq. (|43[) were calculated from Eq. ([2]). 
We assumed that the atoms are specularly reflected at the 
surface of the dust grain (Si — 1). As the initial conditions 
for a dust particle with R — 2 fim, mass density g — 1 
g/cm 3 , and Q' pr — 1, we used a in — 60 AU, e in = 0.2, uj in 



= 120°, Clin — 30°, and ii„ — 20°. The initial true anomaly 
of the dust particle was f in — 180° for Eq. (143)1 . The re- 
sults are depicted in Fig. [5] The solid lines are used for the 
solution of Eq. (|43[l and the dashed lines are used for the 
solution of Eqs. (|44 || - (j48)) . The black lines are used for the 
variable drag coefficients and the grey lines are used for con- 
stant drag coefficients. The solution of Eqs. (|44[l - (|48p with 
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Figure 6. Orbital evolution under the action of the PR effect, radial solar wind, and interstellar gas flow, obtained from numerical 
solution of Eq. H43H (black line) and from numerical solution of the system of differential equations Eqs. I14l| - |48l l (grey line). We used 
interstellar gas with three different temperatures T\ = 500 K (solid line), T\ = 5000 K (dashed line) and T\ = 50000 K (dotted line). 



constant drag coefficients can be obtained by putting fc, ; = 
in Eqs. (|4"4")) - ((4"g)) . Fig. [S] shows that the solution of the 
equation of motion (Eq. I43|) is in good accordance with the 
solution of the system of equations constituted by the sec- 
ular time derivatives of the Keplerian orbital elements Eqs. 
(|44[) - (|48[) . for both variable and constant drag coefficients. 
The semimajor axis decreases faster for the constant drag 



coefficients. This is in accordance with the properties of the 
function gi (see the discussion after Eq. [2D and Eq. l22p . The 
numerical solutions depicted in Fig. [5] represent the cases 
for which Eq. ([8| holds. In these cases, the influence of the 
variable drag coefficient in the acceleration caused by the 
interstellar gas flow on the orbital evolution of the dust par- 
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tide is not large and in some cases can be neglected (as we 
can see in Fig. [5]). 



4.3 Validity of the linear approximation at 
various Mach numbers 

Fig. [6] compares solutions of Eq. (|43[) (black line) with so- 
lutions of Eqs. (|44)| - (j48)) (grey line). The variability of the 
drag coefficient in Eq. (I43|) was given by Eq. (J2J. We used 
an artificial interstellar gas flow which consists only of neu- 
tral hydrogen atoms with concentration n\ =0.1 cm -3 . The 
hydrogen gas velocity vector with respect to the Sun was vf 
= (10 km s , 25 km s , 5 km s _1 ). In order to visualise 
the influence of the molecular speed ratio (Mach number) on 
the orbital evolutions, we used three different temperatures: 
Ti = 500 K (solid line), 2\ = 5000 K (dashed line) and 7\ 
= 50000 K (dotted line). These parameters correspond to 
Mach numbers (the first equation in Eqs. [TTj) soi = 9.5, Sol 
= 3.0, and soi = 1.0. Condition soi <C 1 does not hold for 
none of these values. Therefore, derivation of Eqs. (|44[) - (|48[) 
using the acceleration caused by the interstellar gas flow de- 
scribed by Eq. {TJ is correct (condition for validity of Eq. [6l 
soi <C 1, is not fulfilled). As the initial conditions for a dust 
particle with R — 2 /im, g — 1 g/cm 3 and Q' pr — 1, we used 
a in = 60 AU, dn = 0.2, <jJi n — 120°, Q. in = 30°, and i in = 
20°. The initial true anomaly of the dust particle was fi n 
- 180° for Eq. |g3J. Fig. [5] shows that the orbital evolution 
under the action of the PR effect, radial solar wind, and in- 
terstellar gas flow, is well described by the solution of Eqs. 
(|44p - (|48|) also for the various values of Mach numbers. The 
fact that evolution of the dust particle under the action of 
an interstellar gas flow with a larger temperature is faster 
is caused by proportionality of the secular time derivatives 
to cqi. cqi is larger for an interstellar gas flow with a larger 
temperature (see the first and the second equation in Eqs. 
(fTTj) and Fig. Q] or Appendix [B]). 



5 CONCLUSION 

We have investigated the orbital evolution of a spherical dust 
grain under the action of an interstellar gas flow. The ac- 
celeration of the dust particle caused by the interstellar gas 
flow depends on a drag coefficient which is a well determined 
function of the relative spee d of the dust particle with re- 
spect to the interstellar gas (jBaines. Williams &: Asebiomd 
1 19651 ). We assumed that the acceleration caused by the in- 
terstellar gas flow is small compared to the gravitation of a 
central object, that the speed of the dust particle is small in 
comparison with the speed of the interstellar gas flow and 
that molecular speed ratios of the interstellar gas compo- 
nents are not close to zero. Under these assumptions, we 
derived the secular time derivatives of all Keplerian orbital 
elements of the dust particle under the action of the acceler- 
ation caused by the interstellar gas flow, with linear variabil- 
ity of the drag coefficient taken into account, for arbitrary 
orientations of the orbit. 



If the variability of the drag coefficient is taken into 
consideration in the acceleration, then the secular decrease 
of the semimajor axis is slower. The secular decrease of the 
semimajor axis is slowest for orbit orientations characterised 
by the perpendicularity of the orbital plane to the interstel- 
lar gas velocity vector. The negative secular time derivative 
of the semimajor axis is in this case independent of the ec- 
centricity of the orbit. The secular decrease of the semima- 
jor axis is for a given orbit fastest in the planar case (when 
the interstellar gas velocity vector lies in the orbital plane) 
with the interstellar gas velocity vector parallel to the line 
of apsides. For such orbits with various eccentricities is the 
secular decrease of the semimajor axis fastest for a orbit 
with largest eccentricity. 

Regarding the secular evolutions of the eccentricity, the 
argument of perihelion, the longitude of the ascending node, 
and the inclination, we found that the variability of the drag 
coefficient has a tendency to compensate the influence of the 
terms multiplied by a, see Eqs. H14p - (|17|) . The terms multi- 
plied by a originate from the dependence of the acceleration 
caused by the interstellar gas flow on the velocity of the dust 
particle with respect to the central object. 

If we consider only the influence of the interstellar gas 
flow on the orbit of the dust particle, then the product of 
the secular eccentricity and the magnitude of the radial com- 
ponent of vf measured in the perihelion is, approximately 
constant during the orbital evolution. A simple approxima- 
tive relation also holds between the secular eccentricity and 
the magnitude of the normal component of vf measured at 
perihelion. 

In the special case when the interstellar gas flow veloc- 
ity lies in the orbital plane of the particle and the particle 
is under the action of the PR effect, the radial solar wind, 
and an interstellar gas flow, the orbit approaches the posi- 
tion with maximal value of the magnitude of the transversal 
component of vf measured at perihelion. 

We found, by numerically integrating the equation of 
motion with a variable drag coefficient, that the linear ap- 
proximation of the dependence of the drag coefficient on the 
relative speed of the dust particle with respect to the in- 
terstellar gas is usable for practically arbitrary (no close to 
zero) values of the molecular speed ratios (Mach numbers), 
if the interstellar gas flow speed is much larger than the 
speed of the dust particle. 



APPENDIX A: DERIVATION OF THE 
SECULAR TIME DERIVATIVES OF THE 
KEPLERIAN ORBITAL ELEMENTS 

We want to find the secular time derivatives of the Keplerian 
orbital elements (a, semimajor axis; e, eccentricity; ui, argu- 
ment of perihelion; il, longitude of the ascending node; i, 
inclination). We will assume that the acceleration caused by 
the interstellar gas flow can be used as a perturbation to the 
central acceleration caused by the solar gravity. We use the 
Gaussian perturbation equations of celestial mechanics (cf, 
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e.g., iMurrav fc Dermottl (|l999h , iDanbvl (|l988h l. Therefore 
we need to determine the radial, transversal, and normal 
components of the acceleration given by the second term 
in Eq. (|12[1 . The orthogonal radial, transversal, and normal 
unit vectors of th e dust particle in a Keplerian orbit are (cf., 
e.g.. |Pastorll2009h 

en — (cosfi cos(/ + ui) — sinfi sin(/ + uj) cosi , 
sinf2 cos(/ + uj) + cos Q sin(/ + ui) cosi , 
sin(/ + ui) sini) , (Al) 

eV = (— cosO sin(/ + ui) — sinfi cos(/ + uj) cosi, 
— sin fl sin(/ + uj) + cos Q cos(/ + uj) cos i , 
cos(/ + u;) sini) , (A2) 

ejv = (sinfi sini, — cos Q sini, cosi) , (A3) 

where / is the true anomaly. The velocity of the particle in 
an elliptical orbit can be calculated from 



df d 
~dt ~ di 
esin/ 



(re R ) 
df 



df 



e R + r e T 
1 + e cos f dt dt 



where 



1 + e cos / 



(A4) 



(A5) 



and p — a(l — e 2 ). In this calculation, Kepler's Second Law, 
df /dt = y/Jip/r 2 , must be used. Now, we can easily verify 
that 



(v — vf ) • en = vf a e sin / — vf ■ Sr 

= vf a e sin f — A , 

(y — vf) ■ eV = Vf a (1 + e cos /) — vf ■ St 

— vf o (1 + e cos f) — B , 

(p — Vf) • Sv — — vf ■ e~N = — C , 

where 

a — . 

vf 

Using the notation defined in Eqs. (|A6|) - (|A8|) and Eq. (|A4|I . 
we can write 



(A6) 

(A7) 
(A8) 



(A9) 



v ■ vf — cr vf [B + e( A sin f + B cos /)] 



(A10) 



If we denote the components of the interstellar gas flow ve- 
locity vector in the stationary Cartesian frame associated 
with the Sun as vf — (vfx, vfy, Vfz), then we obtain 



A sin f + B cos / 



Hence, 

V ■ vf — cr vf (B + el) 



(— cosf2 sin a; 

— sinfi cosoj cosi) vfx 

+ (— sinfi sinuj 

+ COS f2 COS UJ COS i) Wfy 

+ cosuj sini dfz = I . (All) 

(A12) 



For radial (an), transversal (ax), and normal (oat) compo- 
nents of the perturbation acceleration, we then obtain from 
the second term in Eq. fl2]l. Eqs. (fA6 |l - ((A"8|> . and Eq. (|A12|) . 
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B f h 

+ a — 1 H vf 
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(A15) 



Now we can use the Gaussian perturbation equations of ce- 
lestial mechanics to compute the time derivatives of the or- 
bital elements. The time average of any quantity g during 
one orbital period T can be computed using 

-i 



(<?) 



g dt = 



W< 



2 n a 3 / 2 



9\i) df 



dt 



y/U 



2 n a 3 / 2 



2 7T a 2 VI -e 2 



2tt 



df 



g r df 



(A16) 



where the Kepler's Second Law, y / 7ip — r 2 df/dt, and Ke- 
pler's Third Law, 47r 2 a 3 = /iT 2 , were used. This procedure 
is used in order to derive Eqs. (|13[) - (|17p . 



APPENDIX B: BEHAVIOUR OF g t 

We define 

l{s) 



gi{s) = 



where 



CDi{s) 
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Us) = — = r I e 



-( l -^+^)"'- t( M (132 ' 



In order to find the behaviour of l(s), we can write 
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dl 2 (s) 
ds 



= 4s 

d_ 
ds 



2s _ s 2 , 
= e + erf(s 



2s _ s 2 N 

— =e + erf(s, 



(B4) 



4s" _ s 2 
-;=e >0 



(B5) 



Since dh (s)/ds ^ 0, h (s) is an increasing function of s for s G 
(0, oo). The value of ^(0) = 0. Therefore h(s) is positive for 
s G (0, oo). If h(s) is positive, then dh(s)/ds > 0. Therefore 
h(s) is an increasing function of s. The value of Zi(0) = 0. 
Thus, h(s) is positive for s G (0, oo). If 'i(s) is positive, 
then dl(s)/ds > 0. Because dl(s)/ds > 0, the function l(s) 
is an increasing function of s for s G (0, oo). lim s ^o l(s) — 
and lim s _ i . 00 l(s) = 1. 

Now, we find the behaviour of coi{s). We can write 



dc Di (s) 

ds 



dc D ii(s) 
ds 



2s _ 3 2 
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(B7) 



Since dcDn(s)/ds ^ 0, cdh(s) is a decreasing function of 
s for s G (0, oo). The value of cdh(0) = 0. Therefore 
Cdh(s) is negative for s G (0, oo). If Cdh(s) is negative, then 
dcDi(s) I 'ds < 0. Because dcoi(s)/ds < 0, the function CDi(s) 
is a decreasing function of s for s G (0, oo). lim s ^o coi(s) = 
oo and lim^oo c D i(s) — 1. 

l(s) is an increasing function of s and cm is a decreasing 
function of s for s G (0, oo). Both l(s) and cd;(s) are positive. 
Therefore, the function gt(s) = l(s)/cDi(s) is an increasing 
function of s for s G (0, oo). 



APPENDIX C: BEHAVIOUR OF h 

We have 



l-VT 



He) = 

In order to find the behaviour of h(e), we can write 

dh(e) _ 2 - e 2 - 2^1 - e 2 



(/e 



3 ^T 



ae oe V 



(Cl) 



(C2) 



(C3) 



Because dhi(e)/de ^ 0, fti(e) is an increasing function of 
the eccentricity. The value of fti(0) is 0. Therefore, hi(e) is 
positive for e G (0, 1]. If fti(e) is positive, then dh(e)/de > 
0. Because dh(e)/de > 0, the function h(e) is an increasing 
function of the eccentricity for e G (0, 1]. 



ACKNOWLEDGMENTS 

I want to thank Francesco Marzari for his useful comments 
and suggestions. 



REFERENCES 

Alouani-Bibi F., Opher M., Alexashov D., Izmodenov V., 

Toth G., 2011, ApJ, 734, 45 
Bahcall J., 2002, Phys. Rev. C, 65, 025801 
Baines M. J., Williams LP., Asebiomo A. S., 1965, MN- 

RAS, 130, 63 
Belyaev M., Rafikov R., 2010, ApJ, 723, 1718 
Burns J. A., Lamy P. L., Soter S., 1979, Icarus, 40, 1 
Danby J. M. A., 1988, Fundamentals of Celestial Mechan- 
ics, 2nd edn. Willmann-Bell, Richmond, VA 
Debes J. H., Weinberger A. J., Kuchner M. J., 2009, ApJ, 

702, 318 
Frisch P. C. et al., 2009, Space Sci. Rev., 146, 235 
Gustafson B. A. S., 1994, Annu. Rev. Earth Planet. Sci., 

22, 553 
Hines D. C. et al., 2007, ApJ, 671, L165 
Lallement R., Quemerais E., Bertaux J.L., Ferron S., 

Koutroumpa D., Pellinen R., 2005, Science, 307, 1447 
Klacka J., 2004, Celest. Mech. and Dynam. Astron., 89, 1 
Klacka J., Komar L., Pastor P., Petrzala J., 2009, in Jo- 
hannson H. E., ed., Handbook on Solar Wind: Effects, Dy- 
namics and Interactions. NOVA Science Publishers, New 
York, p. 227 
Klacka J., Petrzala J., Pastor P., Komar L., 2012a, MN- 

RAS, 421, 943 
Klacka J., Petrzala J., Pastor P., Komar L., 2012b, Icarus, 

submitted 
Maness H. L. et al., 2009, ApJ, 707, 1098 
Marzari F., 2012, MNRAS, 421, 3431 
Marzari F., Thebault P., 2011, MNRAS, 416, 1890 
Murray C. D., Dermott S. F., 1999, Solar System Dynam- 
ics. Cambridge University Press, New York, NY 
Pastor P., 2009, preprint (astro-ph/0907.4005) 
Pastor P., 2012, Celest. Mech. Dyn. Astron., 112, 23 
Pastor P., Klacka J., Komar L., 2011, MNRAS, 415, 2637 
Poynting J. H., 1903, Philos. T. R. Soc. Lond., 202, 525 
Robertson H. P., 1937, MNRAS, 97, 423 
Scherer K., 2000, J. Geophys. Res., 105, A5, 10329 
Whipple F. L., 1955, ApJ, 121, 750 
Wyatt S. P., Whipple F. L., 1950, ApJ, 111, 134 



